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ABSTRACT 

The atmospheric composition and vertical structure of the super-Earth GJ1214b has been a subject of debate since its discovery in 
2009. Recent studies have indicated that high-altitude clouds might mask the lower layers. However, some data points that were 
gathered at different times and facilities do not fit this picture, probably because of a combination of stellar activity and systematic 
errors. We observed two transits of GJ1214b with the Large Binocular Camera, the dual-channel camera at the Large Binocular 
Telescope. For the first time, we simultaneously measured the relative planetary radius k = R p /R * at blue and red optical wavelengths 
(B + R), thus constraining the Rayleigh scattering on GJ 1214b after correcting for stellar activity effects. To the same purpose, a long¬ 
term photometric follow-up of the host star was carried out with WiFSIP at STELLA, revealing a rotational period that is significantly 
longer than previously reported. Our new unbiased estimates of k yield a flat transmission spectrum extending to shorter wavelengths, 
thus confirming the cloudy atmosphere scenario for GJ 1214b. 
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1. Introduction 


The observational study of exoplanetary atmospheres is a very 
young field of research that started in 2002 with th e space-base d 
detec tion of sodium around a gaseous giant (ICharbonneau et aD 
20021) and was followed a few years later by the first such de¬ 
tection from the ground dRedfield et al.1120081) . In this context, 
transiting exoplanets offer a great opportunity to characterize 
their atmospheres. During a transit the amount of stellar light 
absorbed or scattered by the planetary limb (and therefore the 
“effective” planetary radius) is wavelength-dependent and can 
be interpreted through models of different chemic al composition 
and physical structure dSeager & Deming1l2010l) . While sharp 
atomic absorption lines such as neutral Na and K are investi¬ 
gated by means of high-resolution spectroscopy, low-resolution 
spectra and broad-band photometry can be exploited to search 
for wide molecular bands in the near-infrared (NIR), or for 
proces ses affecting the continuum such as Rayleigh scattering 
dSing et aTl 1201 ll) . The latter is more efficient in the ultravio- 


* Based on data acquired using the Large Binocular Telescope (LBT). 
The LBT is an international collaboration among institutions in the 
United States, Italy, and Germany. LBT Corporation partners are the 
University of Arizona on behalf of the Arizona university system; Isti- 
tuto Nazionale di Astrofisica, Italy; LBT Beteiligungsgesellschaft, Ger¬ 
many, representing the Max-Planck Society, the Astrophysical Institute 
Potsdam, and Heidelberg University; the Ohio State University; and 
the Research Corporation, on behalf of the University of Notre Dame, 
University of Minnesota and University of Virginia. Partly based on 
STELLA WiFSIP data (Strassmeier et al. 2004). 

** email address: valerio.nascimbeni@unipd.it 


let (UV) and blue optical regions due to its A -4 dependence. 
The positive or null detection of Rayleigh scattering can greatly 
help in distinguishing between families of atmospheric mod¬ 
els that_otherwise would be strongly degenerate, as shown by 
iHowe & Burrowsl(i2012l) . With the advent of more sophisticated 
instruments, astronomers are focusing on characterizing smaller 
and smaller targets, now reaching a class of planets with masses 
between 2 and 10 M ffi . These planets, dubbed super-Earths, lack 
any analog within our solar system, and their possible compo¬ 
sition is still a matter of debate. The missing parameters in¬ 
clude the mean molecular weight // of their atmospheres, a key 
qu antit y that als o co nstrains the inner structure of these planets 
(IValencia et al .1120131) . 

The first di scovered member of the super-Earth class 
was GJ1214b (ICharbonneau et akl 120091) . a 2.8 /?@, 6.3 
Mc ; low-dens ity p lanet hosted by a M4.5V red dwarf 
dAnglada-Escude et al.ll2013 ). The atmosphere of GJ1214b has 
been investigated with differ ent techniques and inst ruments by 
many authors (among them. Desert et akll201 ll and iBean et alJ 
120111) , sometimes presenting contradictory results. The lack of 
consistency can be explained with a combination of 1) an in- 
comp l ete treatm ent of correlated noise, also called red noise 
dPont et al.ll2007t) : 2) i nhomogeneous ass umptions (for example, 
on orbital inclination: Southworthll2008 :): and 3) a high level of 
stellar activity of the host GJ1214, which has been known for 
a long time to be a spotted and flaring star dBerta et al.ll201 ll) . 
More recent works agree on detecting a flat and featureless 
transmission spectrum, which rules out a cloudl ess, low-/.; atmo - 
sphere as for a solar-composition atmosphere dBerta et al.ll2012b 
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iFraine et al.l2013Llkreidberg et al.l2014lf . Flat spectra can be due 
both to high-altitude clouds masking the lower layers, or to a 
high-// atmosphere consisting, for instance, of methane or water 
vapor. The two model families are very difficult to distinguish by 
standard techniques, although recently there is some consensus 
about models based on cloud s dMorlev et al .120 1 BllBarstow et al.1 
1201 lllKreidherg et al.l20l4) . Most of these works were based on 
NIR spectrophotometry. 

Probing Rayleigh scattering could be useful to constrain p in 
a complementary way, a s was done in our previou s work about 
the hot Uranus GJ3470b dNascimbeni et al.112013hi) . At the same 
time, it would distinguish between whole classes of models that 
are based on mixtures o f clouds and hazes, as seen for instance in 
Fig. 21 from lHowe & Burrows ; d2012i) . The best optical spectral 
region in which to carry out such a search is that encompassing 
the Johnson U and B bands, where the Rayleigh signal is at its 
strongest. Unfortunately, the blue flux of GJ1214, an extremely 
red star with B = 16.4, is too weak to allow high-precision pho¬ 
tometry with medium-sized telescop es. The only previous works 
that investigated t his spectral region dde Mooii et al.l20T2l 120 1 3b 
iNarita et~afll2013l) delivered relatively low-S/N light curves and 
large error bars on the transmission spectrum. On the other 
hand, all those measurements were not simultaneous, which 
gave rise to some doubts that activity effects (which are both 
wavelength- and time-dependent) might have biased the results. 
As our knowledge on stellar activity progresses, it is becom¬ 
ing clearer and clearer that simultaneous multiwavelength tran¬ 
sit observations, to be carried out with differential photometry or 
spectrophotometry, are much more informative than sparse sam¬ 
plings bec ause at least the time-dependent systemat ic component 
is zeroed dBallerini et al.ll20l2l lOshagh et al .112013l) . 

In this work, we present for the first time a simultaneous 
high-S/N measurement of the relative radius k — Rp/R+ of 
GJ1214b in the B and R bands, obtained through Large Binoc¬ 
ular Camera (LBC) dual-band photometry. Stellar activity and 
rotation was monitored during the same season by a long-term 
photometric follow-up with WiFSIP at STELLA. This paper is 
organized as follows. In Sect. 2 we describe how the observa¬ 
tions were performed and how the light curves were extracted. 
In Sect. 3 we estimate the period and amplitude of the rotational 
modulation of GJ1214b, exploiting this information to correct 
our planetary radii derived from the LBC for activity effects. In 
Sect. 4, we finally compare our results with the literature by con¬ 
structing and interpreting a revised transmission spectrum. 

2. Observations and data reduction 

2.1. LBC dual-channel photometry 

We observed two complete transits of GJ1214b during the nights 
of March 29 and May 17, 2012 with the LBC camera mounted at 
the double 8.4m Large Binocular Telescope (LBT). We refer to 
these two events as T1 and T2 hereafter. The LBC is a wide-held, 
prime-focus imager consisting of two independent channels that 
are op timized for the optica l UV/blue and red region, respec¬ 
tively dGiallongo et al.l 12008 ). In particular, the extremely high 
blue efficiency of this instrument, combined with its simultane¬ 
ous dual-channel capability, makes it unrivaled among other 8- 
10m class facilities. The sky was clear and photometric on both 
nights. We mounted a Bessel B and Bessel R filter on the blue 
and red channel, respectively. Exposure times were set to 30 s 
on both channels during the first transit; during the second one, 
the exposure time on the blue camera (where the flux of GJ1214 
is much fainter) was increased to 60 s to reduce overheads and 


increase the duty cycle. Hard defocusing was applied to both 
cameras (up to 2" and 6" FWHM on the B and R images, re¬ 
spectively) to avoid saturation and minimize pixel-to-pixel sys¬ 
tematic errors. Only a 5' x 7' CCD window from a single chip of 
the LBC mosaics was read out to improve speed. Nevertheless, 
the imaged field of view (FOV) included about fifteen reference 
stars with a S/N of the same order of our target. 

The raw images were corrected for bias and sky flat-field 
using standard procedures. Stars drifted by about 50 pixels 
throughout the series. This caused the target to cross a CCD 
row that is larger than the average pixel size; this is a re- 
sult of a well-known m anufacturing defect first described by 
[Anderson & Kind (119991) . This defect required a special back¬ 
ward correction that was applied as prescribed by the same pa¬ 
per. The light curves of GJ1214b were then extracted by run¬ 
ning the STARSKY pipeline, a code optimized to deliver high- 
precision differential photometry over defocused images, origi- 
nally developed for the TASTE project dNascimbeni et al.ll201 ll 
l2013al) . We refer to our previous work on GJ3470b for fur- 
ther details on how STAR SKY has been fine-tuned on the LBC 
dNascimbeni et al.ll2013b ). The four final light curves, hereafter 
referred to as T1 /B, T1 /R, T2 /B, and T2 /R, are shown in the 
left panel of Fig. Q] Their reduced ^ 1 2 , evaluatecQ on the off- 
transit part ranges between^ 2 = 1.1 and 1.5, meaning that our 
noise budget is dominated by white-noise sources. The dominant 
sources of noise are expected to be mostly photon noise in the 
B light curves and photon noise plus scintillation in the R light 
curves dHowellll2006l) . A closer inspection of the effect of corre¬ 
lated noise within the overall noise budget is presented in Sect. 
13.31 The average photometric error of the T1 /R series is only 0.5 
mmag over a 30 s time scale, which makes it one of the most ac¬ 
curate light curves of GJ1214b ever taken from a ground-based 
facility. 

2.2. STELLA photometric follow-up 

We monitored the host star GJ1214 in the observing seasons 
2012 and 2013 to use its photometric variability as an indi¬ 
cator for the spot coverage on the visible hemisphere at the 
time of the transit observations. The program was performed 
with the 1.2m robotic twin telescope STELLA, located in the 
Observatorio del Teide at Tenerife, Canary Islands, and oper- 
ated by the Leibniz-In stitute for Astrophysics Potsdam (AIP; 
IStrassmeier et al.ll2004l) . We used the wide-field imager WiF¬ 
SIP, which consists of a 4kx4k back-illuminated CCD with a 
plate scale of 0.322"/pixel and four read-o ut amp lifiers . It of- 
fers a field of view of 22' X 22' on the sky (iWeber et al .11201 2 !) . 
Between March 21 and October 10, 2012 we obtained data on 
110 nights; and between March 24 and October 4, 2013 we ob¬ 
served on 81 nights. In 2012, we scheduled STELLA to observe 
GJ1214 with blocks of five exposures through a Johnson V fil¬ 
ter and an exposure time of 150 s, followed by five exposures 
through a Cousins I filter with 60 s of integration time. To moni¬ 
tor the host star activity at wavelengths more sensitive to starspot 
modulation, in 2013 we switched to filters Johnson B (three ex¬ 
posures per block, 300 s of exposure time each) and Johnson V 
(three exposures, 150 s). In total, we gathered 197 images in B, 
955 images in V, and 739 in I. 

Bias and flat-field correction was applied by the STELLA 
data reduction pipeline (for details see Granzer et al. 2015, in 
prep.). The I -band data showed fringe patterns whose shape var- 


1 We define^ 2 as the^ 2 divided by the number of degrees of freedom, 

i.e., the number of data points minus the number of free parameters. 
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Fig. 1. Left panel: light curves of GJ1214b gathered through the blue channel of LBC in the Bessel B band (blue points), and through the red 
channel in Bessel R (red points). Data points are plotted with the original sampling cadence. The black line corresponds to the best-fit model 
adopting the parameters listed in Table[2] Right panel: same as above, but showing the residuals from the best-fit model. The bump-like feature at 
phase 0.070 < <f> < 0.092 appeared during the second transit and is interpreted as a starspot crossing; it is marked with an arrow in both channels. 


ied mainly with the CCD temperature. In addition to changes in 
the sky OH emission lines and the sky continuum background, 
a second-order effect in the fringe amplitude is due to a change 
of opacity in the CCD coating when the detector temperature 
changed due to low vacuum, for example. We split all our I - 
band measurements into bins of different CCD temperatures and 
created a master fringe map for each bin as the average of the in¬ 
dividual images after they were filtered from point-like sources 
and smoothed. The master fringe map of the proper temperature 
was then scaled to match the fringe amplitude of each individual 
science frame and subtracted from it. This resulted in a suppres¬ 
sion of the fringe features to about 20% of their original value. 

Aperture photometry was performed with Source ExtractoiQ 
(iBertin & Arnoutsll 19961) . We measured stellar fluxes with the 
Source Extractor parameter MAG_AUT0, which computes an el¬ 
liptical aperture according to the second-order moments of the 
object’s light distribution. This option provides the flexibility to 
account for very different shapes of the PSF caused by differ¬ 
ent observing conditions over the length of our follow-up pro¬ 
gram. We applied stringent selection criteria by rejecting data 
points above a certain threshold in airmass {X — 1.7), sky back¬ 
ground (2 000 ADUs), object FWHM and elongation. Further¬ 
more, we rejected images where the target falls onto the edge 
between different read-out amplifiers. The final sample consisted 
of 128/436/425 images in B/V/I, respectively. To perform differ¬ 
ential photometry, we constructed an artificial comparison star 
by summing the flux of several stable reference stars. We tested 
several selection criteria and found that a combination of four 
nearby stars within 3' from GJ1214, imaged on the same read¬ 
out amplifier, minimized the scatter in the final light curve. The 
same comparison stars were employed for all three data sets B, 
V, and I. 


2 http://www.astromatic.net/software/sextractor 


The final light curves are shown in the two upper panels of 
Fig. [2] in their original cadence (gray dots). Fight curves were 
also binned by averaging all three or five points within a single 
block, to obtain measurements that are much more robust against 
outliers (green, red, and blue circles in Fig. [2] for the V, /, and 
B data, respectively). From here on, our analysis refers to the 
binned light curves. 

3. Data analysis 

3.1. Starspot crossing event 

After a few preliminary fits of the FBC data, it became obvi¬ 
ous that a bump-like feature (marked with an arrow in the bot¬ 
tom right panel of Fig.[T} is clearly detectable in the residuals of 
both channels, just before the egress of the second transit (T2 /B, 
T2 /R). We fitted each residual light curve (i.e., after the transit 
model had been subtracted) with a Gaussian function by leaving 
five free parameters: its height over the transit model 6, its width 
cr, its center to, and the coefficients Ao, A] of a linear baseline 
Ao + A[t. The best-fit results, tabulated on Table |T| and plotted 
in the left panel of Fig. [3] were obtained through a nonlinear 
least-squares algorithm, while the uncertainties were estimated 
though 10 000 iterations of a Monte Carlo bootstrap. 

Within the error bars, the duration and phase of the bump is 
consistent in both light curves, while the maximum amplitude 
is about 40% larger in the B band. This feature can be easily 
interpreted as the crossing of an active region on the stellar sur¬ 
face. Indeed, the contrast between the unperturbed photosphere 
and a colder starspot is higher at shorter wavelengths, giving rise 
to this typical color signature. We emphasize that the two FBC 
cameras are fully independent, and so are the optical paths of 
the two FBT telescopes when observing at the prime foci. Fur¬ 
thermore, no similar feature has been detected in either chan¬ 
nel among field stars with a red color and similar magnitude to 
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Fig. 2. Upper panel: STELLA/WiFSIP differential light curve of GJ1214b, gathered through a Bessel V filter during the 2012 and 2013 observing 
seasons. Unbinned points are shown in gray, daily averages are plotted with green points. Epochs corresponding to the two transits observed by 
LBC (T1 and T2) are marked with vertical lines. Middle panel: same as above, but through a Bessel 1 filter in 2012 (red points) and Bessel B in 
2013 (blue points). Lower panels: GLS periodograms of the light curves, plotted with matching colors. The V periodogram has been evaluated on 
both the whole 2012/2013 data set (solid line) and only on the 2012 season (dashed line). The best-fit period, or the lower limit to it, is reported 
for each periodogram. 


Table 1 . Starspot parameters as fitted by a Gaussian model on the 
residual light curves of the second transit (T2/B and T2/R light curves, 
Fig.[3}. 



72/B 

72/R 

unit 

8 

4.9 ±0.8 

3.3 ±0.3 

mmag 

to 

8.07 + 0.16 

8.12 + 0.09 

orbital phase %o 


44.71 + 0.36 

44.83 + 0.20 

minutes from T\ 


0.848 ± 0.007 

0.851 +0.004 

transit phase 

cr 

1.07 + 0.19 

1.10 + 0.11 

orbital phase %c 


2.44 + 0.43 

2.50 ± 0.24 

minutes 

FWHM 

5.7+ 1.0 

5.9 ±0.6 

minutes 


Notes. The transit phase is cal culated by assuming the total duration 
of the eclipse (Tt-Ts. fkippinglflOl O). in our case 52.4 minutes. The 
FWHM is calculated from the Gaussian cr as 2 V21n 2<x. 


GJ1214 (Fig. [3] left panel, lowest three curves). We can there¬ 
fore reject the hypothesis of an instrumental or telluric origin of 
the bump. 

A few authors already re ported tentative detections of 
starspots occulted by GJ1214b (fBerta et al.l 201 it ICarter et~aT1 
1201 It iKundurthv et al.ll201 lb iNarita et al]|2013l)~ ~but this is the 
first firm detection based on multichannel photometry. The size 
of this active region is probably of the same order as GJ1214b: 


the approximate timescale of the crossing is 5.8 + 0.5 min, which 
is compatible with the time the planet takes to cross a point-like 
feature on its host star (6.2 minutes). We caution that since the 
crossing took place very close (~2 min before) to the third con¬ 
tact of the transit, where the angle between the normal to the 
stellar surface and the line of sight is 6 ^ 65°, the amount of 
photosphere covered by the occulted spot could be up to about 
twice the disk area of GJ1214b (~2%), which is a lower limit for 
the overall spot coverage. 

The more obvious effect of a starspot-crossing event on the 
derived transit parameters is to decrease the best-fit value of 
k, as well as pertu rbing the other parameters in unpredictable 
ways (lOshagh et al.1120131) . To avoid this, in the following anal¬ 
ysis we masked the bump from the 72/B and 72/R light curves 
by assigning zero weight to the affected points, at orbital phase 
0.070 < </> < 0.092. 

By assuming that the size of the spot is the same as the pro¬ 
jected planetary disk, we can obtain a crude estimate of the tem¬ 
perature difference AT between the spot and the unspotted pho¬ 
tosphere. If one compares the heights 6 of the B and R bumps 
(Table |T} with the depth of the underlying transit model at the 
central instant of the spot (14.0 mmag under the out-of-transit 
level), it is straightforward to derive the spot color relative to the 
average stellar color: (B-R) snn ,-(B-R)+ - 0.17 mag. By assum¬ 
ing T e ff* - 3250 K (tAnglada-Escude et al.ll2013l) and modeling 
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Fig. 3. Left panel: modeling of the T2 /R (left plots) and T2/S (right plots) light-curve bumps. From the top: residual light curves with the best-fit 
Gaussian model overplotted (black line; Table [TJ. residuals from the Gaussian fit, and light curves of three reference stars with a red color and 
magnitude similar to GJ1214. Within the errors, the bump-like feature is not seen in any other light curve. Right panel: Correlated noise diagrams 
(photometric RMS of the binned time series as a function of bin size) for the GJ1214b LBC residual light curves after the best-fit transit model 
has been subtracted. The lower sequence in the T2/B, T2 /R plots is evaluated after the starspot feature has been masked out. The green lines follow 
a Poisson-like 1 / VV scaling law as extrapolated from the RMS of the unbinned curve (dotted-dashed line), and with a 50% increase as reference 
(dotted line). 


both fluxes as black-body radiation, we obtain A T - 110 K. 
This agrees with previous observational studies tha t mea sured 
AT < 500 K on early-to-mid M dwarf stars (lBerdvuginall2005l : 
iBarnes et al,ll201 ll iJackson & Jeff'riesl 120 1 3l) . Our AT guess is 
of course a lower limit to the true AT, as in principle the spot 
could be smaller than the occulting planetary disk, and as a con¬ 
sequence, T spot could be cooler than our estimate. 

3.2. Transit parameter extraction 

The LBC light curves have to be fitted with a transit model to 
extract the parameters of interest, above all the relative radius 
k = R p /R+ , which is the wavelength-dependent quantity needed 
to construct the tran smissi on spec trum. We used the JKTEBOP 
cod43 version 34 ( ISouthworth et al.ll2004 ) to perform this task, 
adopting a quadratic law to model the stellar limb darkening 
(LD; [Claret 12004). In principle, the free parameters of the fit 
are the central transit time To, the orbital period P, the fractional 
stellar and planetary radii R+/a and R p /a, the orbital inclination 
i, the Ao, A] coefficients of a linear baseline Ao + Aik and the 
linear and quadratic LD terms u\, no. Within this model, k is a 
derived parameter obtained by dividing R p /a by //*/«. 

We did not fit for both LD terms, following a commo n prac¬ 
tice to improve the robustness of the fit (ISouthworthl 120081) . 
Instead, U 2 was fixed on all light curves at its theoretical 
value for the corresponding photometric band by interpolat¬ 
ing the LD tables by iCIaret et~ahl ( 20121) computed for cool 
stars (1500 < T eff < 4800 K). We adopted for GJ1214 
the most recent v alue s of T e s, log# and [M/H] estimated by 
lAnglada-Escude et al.l ( 2013 ). We investigated the effect on our 

3 http://www.astro.keele.ac.uk/~jkt/codes/jktebop.html 


fit of the linear term u\ by setting it as a free parameter, or by 
keeping it fixed at its theoretical value over independent JK¬ 
TEBOP runs. GJ1214b does not show any tra nsit time varia¬ 
tion (TTV) within about 10 s (iHarpsele et al.l l20l3i) . therefore 
we fixed_both Tp and P at their most recent and accurate values 
dKreidberg et alJl2014l) . Furthermore, the orbi tal inclination i is 
well k nown to be strongly correlated with k (ISouthworth et al.l 
12005 ). To avoid biases when comparing our results with the lit¬ 
erature, we fixed the inclination and scaled semi-a xis at their 
best-fit value i = 88°.94, a/R+ — 14.97 found by iBean et all 


eluding Bean et al. 20 

1; Desert et a 

. 2011 

Croll et al. 

2011 

de Mooii et al. 20121 

Murgas et al.l 

20l2 

Teske et al.i 

2013 


Narita et al.ll2013: Fraine et al.ll2013l.ide Mooii et al.ll2oT3 


After these considerations, we had only three or four fitted 
parameters for each light curve: k, Ao, Ai and, when left free. 
Mi. The other para meters were injected into the fit as Gaussian 
priors, as done in iNascimbeni et al.l (l2013bl) : t heir dispersions 
were s et to the published error s for P, To and i bv lKreidberg et al.l 
(1201 4l) and IBean et aD (12010l) . and as the errors propagated on 
u i and «2 by assumi ng the published uncertainties on 7/n, log# 
and [M/H] bv lAnglada-Escude et al.l (12013l) . We chose to derive 
the uncertainties over the best-fit param eters with a Monte Carlo 
bootstrap algorithm (JKTEBOP task 7: ISouthworth et al1l2005l) 
because formal errors are known to be underestimated due to 
correlations between p aramete rs and residual systematic errors 
(the so-called red noise: lPont et al.|[2007t a detailed analysis of it 
follows in the next subsection). The results are plotted in Table 
[2] the four individual and the two combined light curves; in the 
latter case, each pair of transits observed through the same filter 
was fitted simultaneously. 
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We note that when leaving u\ free, the best-fit value of the 
first epoch is consistent with the theoretical prediction. The best- 
fit values of k are also consistent with those obtained by fitting 
the light curves of the second transit in the matching filter. For 
this reason, and because the masking of the bump-like feature 
does not appear to bias the k measurement, the simultaneous fit¬ 
ting of both transits appears justifiable. 


3.3. Correlated noise analysis 

After the transits were properly modeled, we investigated 
whether and to which extent the residual light curve (i.e., af¬ 
ter the transit best-fit model has been subtracted) is affected by 
correlated noise. To this purpose, a commonly used tool is the 
so-called correlated noise diagram (CND), where the input data 
are binned on different timescales and the resulting photometric 
RMS of the binned t i me series is plotted a s a function of bin 
size dPont et al.ll2007t ICarter & Winnll2009l) . We plotted CNDs 
for all our four light curves (T1 /B, 11/R, T2 /B, and T2 /R) in the 
right panel of Fig. [3] For the 12/B and 12/R transits, we sepa¬ 
rately plotted the CNDs for the full light curve (upper sequence) 
and after masking the starspot feature (lower sequence). The 
CNDs have to be compared with a Poissonian 1 / y/N scaling law 
(where N is the number of points within each time bin), which 
is expected to predict the RMS after the averaging process in 
the ideal case when only white noise contributes (dashed-dotted 
green line on Fig. [3]). More realistically, where correlated noise 
is not negligible, we expect to see an increase of the RMS to¬ 
ward the longest timescales. This is clearly visible in our case. 
For the 11/R and 11/B light curves, the RMS increase is always 
within ~20% of the overall budget even at the longest avail¬ 
able timescales (bin size ~20 min, i.e., the maximum allowed 
by our requirement of having a set of at least ten bins to evaluate 
a meaningful RMS). This proves that the effect of red noise on 
the first transit is marginal. On the other hand, the 12/R and 12/B 
light curves show more red noise, which boosts the binned RMS 
by ~50% even if the starspot bump is not taken into account. This 
is also confirmed by the best-fit linear slopes, which are steeper 
for the first transit: <x(Tl/R) oc A -035 , cr(Tl/B) oc A~ a33 , with 
respect to the second one: cr(T2/R) oc A~ a23 , cr(T2/B) oc AT 0 - 22 
(the latter values are calculated after masking the spot). 

Although the T2 light curve is still dominated by photon 
noise at the 10-15 min timescale, the question is whether this 
additional noise is of instrumental, atmospheric, or astrophysical 
origin. A closer inspection of the 12/R vs. 12/B residual curves 
(lower right panel of Fig. 0 reveals some correlation between 
the two channels, but only during the transit. To support this hy¬ 
pothesis, we selected three subsets of data points covering the 
off-transit part, the in-transit (spot excluded), and the spot feature 
itself (Fig. [4j upper panel). Points were then paired by matching 
mid-exposure times that differ by less than 30 s between the two 
channels. By comparing the B and R residuals of the spot and 
in-transit series, a strong positive correlation is found (Spear¬ 
man rank correlation coefficient: g — 0.833 and 0.721, respec¬ 
tively), while the off-transit points are essentially uncorrelated 
(g = 0.023). Interestingly, if we fit a straight line to the first 
two subsets (Fig. [4] lower panels), even the slope of the corre¬ 
lation is similar (0.58 vs. 0.46; i.e., the B amplitude in on aver¬ 
age about twice the R one). As a similar behavior is not seen in 
any reference star, we conclude that most of the measured red 
noise is due to the crossing of additional inhomogeneities on the 
GJ1214 spotted surface, which was caught here in a relatively 


active phase. Additional evidence of such an increase of activity 
is exposed in the next subsection. 

To check whether and to which extent the correlated noise 
content of our data might bias our derived transit parameters 
(Table 0. we performed an additional error analysis based on 
a residual-permutation (RP) or prayer bead algorithm, which is 
based on the cyclic permutation of the residuals (ISouthworthl 
120081 ). This approach preserves the temporal structure of the cor¬ 
related noise and propagates it through the whole light curve, 
leading to more realistic error estimates. We ran the algorithm 
on our combined data sets, setting JKTEBOP with the same 
input parameters described in the previous subsection. The re¬ 
sulting best-fit values are k = 0.1176 ± 0.0009 (B) and k - 
0.1175 ±0.0004 (R), which are fully consistent with the previous 
ones but with a slightly larger error bar. Following a conservative 
choice, we adopted the RP best-fit values to build the transmis¬ 
sion spectrum of GJ1214b. 

3.4. Stellar activity and rotation 

All our three B, V, I STELLA light curves of GJ1214 (Fig. 0 
show clear signs of periodic variability at similar timescales. 
The phase of that signal is coherent among different filters (V 
vs. I in 2012 observing season, B vs. V in 2013), while the am¬ 
plitude is larger at shorter wavelengths. This signature is easily 
interpreted as rotational modulation, meaning that it is the re¬ 
sult of an unevenly spotted photosphere on a rotating star. The 
photometric variability and flaring behavior of GJ1214 is ex¬ 
pected forjj_M4 1 5Vjlw'arfand was already known before our 
study dKundurthv et alJ 1201 ll). but the rotationa l period P mt is 
not firmly constrained so far. iBerta et ahl (1201 ll) tentatively de¬ 
tected a signal with an amplitude of A rot = 3.5 ± 0.7 mmag 
and P rot - 53 days from their MEarth (i + z)-band light curves 
gathered in 2010, but the authors themselves cautioned that their 
other data sets (Mearth 2008, 2009 seasons and FLWO 2010 V- 
band photometry) are not able to c onfirm that_signal. A more 
recent multicolor follow-up bv lNarita et al.l (1201 3l) carried out at 
the MITSuME 50 cm telescope in 2012 found P rot = 44.3+1.2 d, 
with A rot = 21 +4,5.6+0. 8 , and 3.2±0.4 mmag on the g' , R c , and 
I c bands, respectively. In this case, the authors cautioned that the 
g'-band detection is weak because of low S/N and that the lim¬ 
ited baseline of their observations prevented them from probing 
longer periods. More in general, it is possible that the true P TOt 
of a star is an integer multiple (2:1 or 3:1) of the detected photo¬ 
metric period because there ar e multiple groups of ac tive regions 
at different stellar longitudes (IWalkowicz et al.f2013l) . 

We analyzed our B, V, and I light curves, binned 
over each observing block of three or five repeated expo¬ 
sures, through _a_ generalized Lomb-Scargle periodogram (GLS; 
IZechmeister & Kursterl2009l) . The error bar of each binned point 
was estimated as the root mean square (RMS) of the single mea¬ 
surements. For the V series, which spans two observing seasons, 
we analyzed both the whole series and each season separately. 
This is useful because typical active regions are expected to 
evolve over timescales ranging from weeks to months, which 
means that the phase coherence could be broken or attenuated 
when folding very long time series. The resulting periodograms 
are plotted in the lower panel of Fig. 0 and the corresponding 
best-fit GLS solutions in Table 0 We emphasize that the overall 
shape of the light curves and the resulting periodograms are not 
influenced by a particular choice of reference stars, as confirmed 
by further checks. 

All three light curves show long-term (~80 d) periodicities 
detected at high confidence. The V whole-series periodogram 
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Table 2. Summary of the best-fit transit parameters of GJ1214b from LBC photometry. 


transit 

band 

N obs 

cr [mmag] 

x 1 

Xr 

U\ 

k = P p /P* 

k' - k + Ak max 

Tl, 2012/03/29 

B 

185 

1.55 

197.3 

1.11 

0.60 ± 0.07 

0.1158 + 0.0013 

— 

Tl, 2012/03/29 

B 

185 

1.56 

200.2 

1.12 

0.5259 1 ' 

0.1170 + 0.0009 

0.1159 + 0.0009 

Tl, 2012/03/29 

R 

207 

0.57 

282.2 

1.39 

0.42 ± 0.02 

0.1170 + 0.0004 

0.1159 + 0.0004 

T2, 2012/05/17 

B 

64 

1.32 

64.0 

1.05 

0.5259 1 ' 

0.1186 + 0.0010 

— 

T2, 2012/05/17 

R 

123 

0.76 

266.4 

2.22 

0.4326 1 ' 

0.1178 + 0.0004 

— 

combined 

B 

249 

1.50 

273.1 

1.11 

0.5259' f 

0.1176 + 0.0007 

0.1165 + 0.0007 

combined 

R 

330 

0.65 

575.5 

1.76 

0.4326 f 

0.1175 + 0.0003 

0.1164 + 0.0003 

combined, RP 

B 

249 

1.50 

273.1 

1.11 

0.5259 f 

0.1176 + 0.0009 

0.1165 + 0.0009 

combined, RP 

R 

330 

0.65 

575.5 

1.76 

0.4326 f 

0.1175 + 0.0004 

0.1164 + 0.0004 


Notes. The columns list: the ID code and evening date of the transit (or combined when the two transits were fitted together), the photometric 
channel (Bessel B or R), the number of data points analyzed, the RMS <x of the residuals from the best-fit model, the x 2 and reduced x 2 of the 
residuals, the linear LD coefficient ui (marked with t when it is input as a Gaussian prior around its theoretical value from lClareD+jiL| |20j_2f). 
the relative radius k of the best-fit model, and the relative radius k' after the correction for unocculted starspots has been applied iBaherinTet al .1 
1201 2l: see Sect. [53] for details). The reported^ 2 and^ 2 are evaluated from the nominal error bars once the bump-like feature has been removed 
(see Sect. l3.lt . The RP error analysis is based on a residual-permutation algorithm (see Sect. l3.3l for details). 


S T2/B light curve (blue) 

" 5 T2/R light curve (red) 
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-0.02 -0.01 0 0.01 0.02 

orbital phase 



Fig. 4. Statistical correlation between the blue (T2/B) and red (T2 /R) LBC light curve of the second transit. Upper panel: Observed minus calculated 
(O - C) residuals from the best-fit transit model as a function of orbital phase for the T2/B (square points) and the T2/R (round points) time series. 
Data are color-coded according to whether they are labeled as off-transit, in-transit (the spot feature being excluded) and spot only. Only points 
whose mid-exposure time can be matched within 30 s to a frame of the other channel are considered here. Lower panels: Correlation plots between 
residuals from T2/B and T2 /R light curves for the three samples defined above. The best-fit straight line calculated by ordinary least-squares is 
overplotted. The Spearman rank correlation coefficient is tabulated for each subset, along with the blue and red RMS. 


shows an excess of power between 60 and 100 d (with maximum 
power at P ^ 83 d), which turns into a single prominent peak at 
P ^ 69 d when only the 2012 season is taken into account (i.e., 
that with the longest coverage). The formal false-alarm probabil¬ 
ity (FAP) is 10~ 13 and 8-10 15 , respectively. The I series reveals a 


signal at P ^ 79 d with an even higher level of confidence (FAP: 
7 ■ 10~ 16 ). Through both V and I filters, the pseudo-sinusoidal 
shape of the signals and their phase coherence are clearly visible 
in the light curves. A curvature perfectly consistent with the V 
light curve gathered in the same season is noticeable for the B 
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Table 3. Most significant photometric periodicities of GJ1214 detected 
in STELLA/WiFSIP light curves by the GLS periodogram. 


filter 

season 

P (d) 

A (mmag) 

FAP 

Vj 

2012-2013 

83.0 + 0.8 

11.4+1.1 

1 ■ io- 1J 

Vj 

2012 

69.0 + 2.0 

15.1 ± 1.3 

8 ■ 10" 13 

Ic 

2012 

79.6 + 2.5 

12.0+1.2 

o 

1 

Os 

Bj 

2013 

> 70 

> 15 

— 


Notes. The columns list: the bandpass and observing season of the pho¬ 
tometric series, the period P and amplitude A of the GLS best-fit solu¬ 
tion, and the formal false-alarm probability (FAP). The Bj solution is a 
lower limit (see Sect. l3.4l for details). 


series. As the smaller baseline of the B data (~100 d, although 
only ~60 d are densely sampled) prevents us from detecting sig¬ 
nals at P > 70 d, we assumed P =*■ 70 as a lower limit. All 
STELLA data therefore converge toward a rotational period of 
Pun - 80 d, or an integer multiple of it. Interestinglv. lBerta et al.l 
(1201 ll) reported P lot = 81 d as the most significant period in their 
MEarth 2008 data, which matches our results very well, with an 
FAP of 10 4 . On the other hand, the 44-d period reported by 
[Narita etldl (1201 3l) could be interpreted as a 2:1 harmonic of our 
peak, caused by two groups of active regions present at that time. 

The V, I photometric amplitudes measured by our group dur¬ 
ing the 2012 campaign (i.e., when the LBC transits were gath- 
ere d) are typic a lly tw o to t hree times larger th an those reported 
by iBerta et ail (1201 ll) and iNarita et~akl (120131) . This suggests a 
phase of stronger activity undergone by GJ1214 in 2012. This 
hypothesis is also supported by the occulted starspots in our 
T2 transit, and in t h e August 12, 2012 light curve published by 
INarita etldl (1 201 3l) . iBerdvuginal (120051) and IStrassmeierl (120091) 
indicated that for an M4 star the typical temperature contrast 
and filling factor of starspots are AT < 500 K and < 20%, 
respectively. This agrees with our estimate of AT > 110 K 
based on the B - R color of the starspot occulted during the T2 
transit, as shown in Sect. 13.II Using the BT-Settl atmospheric 
model (IAllard et al.ll201 ll) to simulate such an M4 spotted pho¬ 
tosphere under the above assumptions for AT and filling factor, 
we can indeed consistently recover the variability amplitude of 
the STELLA photometric series. 

3.5. Correction for unocculted starspots 

To construct a transmission spectrum of GJ1214b, the key quan¬ 
tity to derive from our LBC transit light curves is the radius ra¬ 
tio k = Rp/R*. This quantity is directly related with the transit 
depth and can be affected by stellar activity in two complemen¬ 
tary ways. The first, most obvious one is by occulted starspots, 
which give rise to bumps in the light curve that force the fit¬ 
ting procedure to underestimate k. Provided that the S/N of the 
timeseries is high enough, these bumps are easily identified and 
masked by assigning them zero weight, as we did on our T2 tran¬ 
sit. 

A second, subtler way to add a bias on k is due to unocculted 
spots, because they decrease the mean surface brightness of the 
visible photosphere, therefore leading to overestimating k. This 
effect is color-dependent (being stronger at shorter wavelengths, 
where the contrast between the unspotted photosphere and the 
cooler spots is increased) and requires a correction that takes into 
account the spot coverage of th e star at the time of observation. 
We followed the approach of iBallerini et al.l (12012l) to compute 
that correction and refer to our previous work on GJ3470b for 


further details dNascimbeni et al.ll2013bl) . We neglected the con¬ 
tribution from faculae, as no typical signatures of facular occul- 
tations by GJ1214b are recorded either in our data or in the liter¬ 
ature light curves. Below, we focus on investigating the impact 
of activity on k during the first LBC epoch (Tl), where we can 
reasonably assume that there are no occulted starspots (case 2 
above). The starspot-crossing event on T2 and the possible pres¬ 
ence of additional minor active regions (see Sect. [331) instead 
complicate the modeling, because a mixed scenario with both 
occulted and unocculted starspots cannot be approximated by 
either of the simplified cases above. 

At the time of Tl transit, the star appeared AV - A/ - 0.023 
mag fainter than the maximum luminosity as it is modeled on 
the U/2012 and 7/2012 light curves by the best-fit sinusoids 
plotted in Fig. [2] Assuming that the net observed stellar flux 
is the sum of the quiet and spotted photosphere, weighted by 
the corresponding coverage factor over the visible hemis phere 
dAfram & Berdvugi~nall2015[ IScandariato & Micelal 12014 ). this 
translates into an estimated sp ot co verage of ~0.020 at the epoch 
of observation. Following lBallerini et al.l (120121) . adopting AT = 
500 K and estimating A = 1 - F sr ,nt/F+ by means of BT-Settl 
atmospheric models (IAllard et al.l 1201 ll) . the maximum correc- 
tiorQ we compute is 


A k mdX (B)/k = -0.0098 

Ak max (R)/k = -0.0097 ( l) ’ l J 

that is, in absolute terms, Ak max (B) = -0.00115 and Ak max (R) - 
-0.00114. 

These corrections can then be applied to the apparent B 
and R radius ratios k(B), k(R) to obtain lower limits for them, 
named k'(B), k'(R) and reported in the last column of Table 
U These corrections might appear rather significant, as they 
change k by about 1%. Nevertheless, the color signature intro¬ 
duced by stellar activity is tiny, because the differential correc¬ 
tion Ak max ( B ) - Ak max (R) < 0.0001 is well within the measure¬ 
ment error bars crfk, It) = 0.0009 and cr(k, R) = 0.0004, respec¬ 
tively. We emphasize that Ak max represents a maximal correction 
and that the true values probably lie somewhere between k and 
k'. 

By comparing k'{B) = 0.1165 ± 0.009 and k\R) = 0.1164 ± 
0.0004, it is straightforward to conclude that the radius of 
GJ1214b is not dependent on wavelength within the error bars, 
even after the effect of unocculted starspots has been accounted 
for. In other words, based on our LBC data alone, where both the 
time- and wavelength-dependent activity effects are corrected 
for, the continuum of the transmission spectrum is flat within 
<1% between It and R. In the next section, we complement 
our measurements with others published in the literature through 
several bandpasses. As for the latter, only the apparent values of 
k are available because most authors of previous studies did not 
apply any correction for stellar activity. For consistency, we are 
forced to adopt our uncorrected k values throughout the com¬ 
parison. However, we demonstrated above that the color effect 
given by the maximal correction (computed on the bluest and 
most affected band of the whole set: LBC B) is negligible with 
respect to the measurement errors, even during a phase of high 
stellar activity (2012). We are therefore confident that such a cor¬ 
rection -even if it were available- would change the mean level 
of the transmission spectrum, but not the slope, leaving the main 
scientific result unchanged. 

4 i.e., when no s tarsp ots are occulted during the transit: case with / occ = 
0 IBallerini et al.li2012t) 


Article number, page 8 of [ 1 1 1 

























































Nascimbeni et al.: An LBT view of the atmosphere of GJ1214b 


It is worth noting that the corrections computed above ne¬ 
glect the contribution from the fraction of spots that are always 
visible during a rotation cycle (because they are located in a po¬ 
lar region) or those that are uniformly distributed over the sur¬ 
face. The latter case is the most likely scenario, as fully con¬ 
vective stars at M* < 0.35 M 0 , such as GJ1214, are expected 
to drive activity in a turbulent dyn amo regime rather than with a 
classic crO dynamo (Barnes et al.l201 ll and references therein). It 
is easy to check what would happen if the filling factor due to the 
component of homogeneously distributed spots is the maximum 
allowed by previous studies for an M4V star (~0.20: lBerdvuginal 
120051) . Models predict that the correction would scale approxi¬ 
mately as a linear function of the filling factor, that is, it would 
be about ten times larger: 

A k max (B)/k = -0.1173 

Ak max (R)/k = -0.1150 {l) 

that is, in absolute terms, Ak max ( B ) = -0.01378 and Ak max ( R ) 

0.01351. We conclude that even after applying such a maximal 
correction Ak max (B) - Ak max (R) - 0.0003, the resulting k'(B) 
and k'(R) would be consistent with each other within the error 
bars, leaving the main scientific result unchanged. 


4. Discussion 

We reconstructed a transmission spectrum of GJ1214b by merg¬ 
ing the two values of k derived in this work (magenta circles 
on Fig. [5} with all the published data points that assumed the 
same prior as we did on orbital inclination and scaled semi-axis 
(centered on / = 88°.94 and a/R* = 14.97, first adopted by 
iBean et af]l2010l) . As explained in Sect. 13.21 this is to avoid a 
systematic bias on k, because k and i (and, to a much lesser ex¬ 
tent, a/R *) are known to be strongly correlated fit parameters. 
Unfortuna tely, this choice ex clu des some very high S/ N data sets 
('includinu lBerta et al.l2012l and lKreidberg et al.l2014l who fixed 
i at 89°.30 and i = 89°.1, respectively), but a full homogeneous 
reanalysis of these data is beyond the scope of this paper. The 
literature data we gathered include the following light curves, 
which resulted in 40 data points after combining measurements 
gathered through the same instrument and setup: 

- one transit from FORS 2 at the VLT ( spectrophotometry 780- 
1000 nm; SP) by IBean et al.l (f20ul) . as later reanalyzed by 
IBean et al.l(l2011l) : 

- three transits from FORS2 at the VLT (610-850 nm SP), 
MMIRS at the Magellan ( H + K SP) and HAW KI at the VLT 
(narrow filter NB2090) bv lBean et alJ (1201 lh : 

- two transits from IRAC at the Spitz er (IRAC 3.6 pm and 
4.6 pm bands) bv lDesert et all (1201 ill : 

- four transits from WIRCam at the CFHT (each light curve 
was gathered by alternatively switching between pair of 
filters chosen among /, K s , and a narrow CFLOn), by 
ICroll et aLld2QTlh : 

- two transits from WFC at the INT (Sloan r and /), one tran¬ 
sit from NOTCam at the NOT (K s ), one transit from LIRIS 
at the WHT (K c ), and one transit simultaneously observed 
through four biters (Sloan g, r, i, z) from GROND at the 
MPI-2.2m, bv lde Mooii et alJ (12012i) : 

- seventeen transits from IRAC at the Spitzer (fourteen at 
4.5 pm and three at 3.6 pm) and seven ground-based transits 
from TRAPPIST (through a wide band encompassing I + z') 
bv lFraine et al.l (120131) : 

- three transits from ACAM at the WHT (Sloan g), FORS2 at 
the VLT (Bessel B) and WFC at the INT (Sloan g) respec¬ 
tively, bv lde Mooii et al.l (l2013h : 


- four transits from Mont4k at the Kuiper-1.55m (two through 
Harris V and two through Harris R) and five transits from 
WiFSIP at STELLA (Sloan g'), bv lTeske et all (120131) : 

- two transits from SuprimeCam at the Subaru and FOCAS at 
the S ubaru (both in Bessel B), and one transit simultane¬ 
ou sly observed throug h J, H, K from SIRIUS at the IRSF, 
bv lNarita et al. ( 2013 ). 


Theoretical models of different c hemical and physical com - 
positions, calculated for GJ1214b by iHowe & Burrows! d2012l) . 
are plotted with solid lines of different colors in Fig. 0 along 
with the mentioned data. We focus on three representative cases: 
1) a low p, cloud-free atmosphere with solar composition, that is, 
dominated by hydrogen and helium (red line); 2) a high p, cloud- 
free atmosphere composed of pure water vapor (H 2 O; blue line); 
3) a featureless (constant) spectrum simulating a high-altitude 
layer of opaque clouds (green line), which mask every spec¬ 
tral feature behind that layer. We integrated the model spectrum 
over the photometric passband of each data point to calculate the 
residuals between model and observation, then we performed a 
statistical analysis. The corresponding^ 2 values per degree of 
freedom are reported for each case in the corresponding diagram 
of Fig. [5] 

The solar-abundance model is clearly ruled out at high confi¬ 
dence, the resu lting y 2 being 7041. A si mil ar conclusion was al - 


dence, the resu iting y being YU4 
ready drawn by|Bean_etalJ(l201L[ 


201 lh and lDesert et al.l (1201 lh . 


debated by ICroll et all ( 201 ll) . and later confirmed by most 


was * 

recent studies. The steep increase of absorption at the blue end of 
this particular model spectrum, which is a signature of Rayleigh 
scattering by molecular hydrogen H 2 , is ruled out at 3<x even by 
the LBC data alone, which makes ours an independent confirma¬ 
tion. 

The water-based model returns a better fit, although with a 
rather high value ofy 2 = 319, while the simpler featureless spec¬ 
trum yields^ 2 = 200. Based on this comparison, we selected the 
latter as the model that best explains the observations. We em- 
ph asize that no other model atmosphere among those computed 
bv lHowe & Burrows! (12012l) returns a smaller^- 2 , while the best- 
fit^- 2 found above implies error bars underestimated by a factor 
of ~2, under the hypothesis of an underlying Gaussian distribu¬ 
tion. This is frequently observed when compar ing results origi¬ 
nating from different studies dSouthworthl201 ll among others, is 
a rich source of such examples). Possible reasons include the use 
of poorly suited numerical algorithms to estimate the uncertain¬ 
ties over the parameters, which cannot be calculated without tak¬ 
ing into account the presence of highly correlated fit variables. 
A second explanation is stellar variability, which can alter most 
transit parameters (including k, the most relevant one to our pur¬ 
poses) as a function of both time and wavelength. For GJ1214b, 
we expect variability to be the dominant source of scatter in its 
optical transmission spectrum, following the discussion made in 
Sect. 13.51 where we estimated that the peak-to-valley relative 
variation of k owing to stellar activity could reach a few percent 
under reasonable assumptions. 

It is more challenging, however, to explain why the apparent 
radius of GJ1214b changes so much when measured at different 
epochs in the NIR, where spot contrast is expected to be lower 
and the corresponding effect on the apparent k weaker. This is 
most strikingly noticeable in the K band, where different au¬ 
thors found radii discrepant by more than 4cr dNarita et al .120131) . 
On one hand, little is known about the detailed magnetic behav¬ 
ior of stars close to the M5V transition where GJ1214 lies, and 
the r eff vs. r spo , relation reported by iBerdvuginal d2005l) stops 
at T e ff = 3300 K. It has been hypothesized that much colder 
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• DeMooij2013 (#6) 
■ DeMooij2012 (#7) 
A Bean2011 (#8) 


Fig. 5. Transmission spectrum of GJ1214b reconstructed by merging results from the present work (magenta circles) with all the published data 
that assumed the s ame prior as we did^on i and a/R * (see the legend on the right). Theoretical models with different compositions calculated 
for GJ1214b by iHowe & BurrowS l l2012i) are plotted with solid lines. Upper panel: low p, cloud-free atmosphere with solar composition (red 
line). Middle panel: high p, cloud-free atmosphere composed by pure water (blue line). Lower panel: featureless (constant) spectrum simulating 
a high-altitude layer of opaque clouds (green line). The corresponding^ 2 values are reported for each case. 


spots could begin to appear at those spectral types dBarnes et al.1 
1201 ll ). That would explain why many surveys of mid- and late-M 
stars measured a photometric jitter th at is sim ilar to that usually 
reported in the optical region dPlavchan et al.li2008l among oth¬ 
ers). 

On the other hand, another plausible interpretation for these 
cases is the presence of residual red n oise, which js known to 
plague NIR data, often in subtle ways (ISwain et al.ll2013l) . 

5. Conclusions 

We demonstrated that the dual-channel capabilities of the LBC 
prime-focus camera and its unrivaled efficiency (especially on 
the blue side of the spectrum) are extremely effective in deliver¬ 
ing high-precision photometric time series of faint, red, and ac¬ 
tive targets such as GJ1214b. Indeed, the Bessel B and R light 
curves presented in this work are the most precise ever pub¬ 
lished in their respective spectral regions. Not less important, 
they are truly simultaneous measurements, which enabled us to 
remove the time-dependent component of stellar activity. We ex¬ 
ploited a smaller, robotic facility (WiFSIP at STELLA) to moni¬ 
tor the photometric variability of GJ1214 over a larger time scale 
to determine the wavelength-dependent contribution. This re¬ 
sulted in a new, robust estimate of the stellar rotational period 
at P rot - 80 d. By comparing the WiLSIP light curves with those 
reported in the literature, we also noted that the activity level 


of GJ1214 was significantly higher than the average through¬ 
out the 2012 observing season. The detection of a large bump 
in the second transit (T2; 2012/05/17), induced by the crossing 
of a starspot, appears to support that finding. The nature of this 
bump is unambiguously confirmed by its B - R signature. The 
latter allowed us to derive a lower limit to the temperature con¬ 
trast between the photosphere and the spot (AT ^ 110 K), again 
highlighting the advantages of multicolor photometry. 

The relative radii of GJ1214b through the B and R filter 
from the LBC data alone, once corrected for the presence of un¬ 
occulted spots, are perfectly consistent with each other within 
~ 1%. This finding suggests a flat transmission spectrum in 
the optical region and rules out at 5cr the tentative detection 
of Rayleigh sca ttering by IT that has been claimed by some 
previous works (Ide Mooii et al.ll2012t iFraine et akll2013l) . After 
combining our points with other optical and NIR measurements, 
we rule out cloud-free atmospheric models dominated by H/He 
and by heavier compositions such as pure water vapor. We con¬ 
clude that a flat transmission spectrum, such as that resulting 
from an atmosphere with a thick cloud layer at high altitude, is 
still the model that best explains the available data. Our conclu¬ 
sion i s co n siste nt with other recent independent works (most no¬ 
tably, j&eidbergjet^OlloIS) that also ruled out cloud-free, high- 
p models in favor of a featureless spectrum. We have to await 
the commissioning of future space-based instruments (such as 
the James Webb Space Telescope) to further investigate the at- 
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mosp here of GJ1214 and confirm the cloudy scenario dBelu et alJ 
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